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Abstract. We study the stability of mixtures of highly screened repulsive charged 
spheres and non-adsorbing ideal polymer chains in a common solvent using free volume 
theory. The effective interaction between charged colloids in an aqueous salt solution 
is described by a screened- Coulomb pair potential, which supplements the pure hard- 
sphere interaction. The ideal polymer chains are treated as spheres that are excluded 
from the colloids by a hard-core interaction, whereas the interaction between two ideal 
chains is set to zero. In addition, we investigate the phase behaviour of charged colloid- 
polymer mixtures in computer simulations, using the two-body (Asakura-Oosawa pair 
potential) approximation to the effective one-component Hamiltonian of the charged 
colloids. Both our results obtained from simulations and from free volume theory 
show similar trends. We find that the screened-Coulomb repulsion counteracts the 
effect of the effective polymer-mediated attraction. For mixtures of small polymers 
and relatively large charged colloidal spheres, the fluid-crystal transition shifts to 
significantly larger polymer concentrations with increasing range of the screened- 
Coulomb repulsion. For relatively large polymers, the effect of the screened-Coulomb 
repulsion is weaker. The resulting fluid- fluid binodal is only slightly shifted towards 
larger polymer concentrations upon increasing the range of the screened-Coulomb 
repulsion. In conclusion, our results show that the miscibility of dispersions containing 
charged colloids and neutral non-adsorbing polymers increases, upon increasing the 
range of the screened-Coulomb repulsion, or upon lowering the salt concentration, 
especially when the polymers are small compared to the colloids. 



1. Introduction 

Adding non-adsorbing polymers to colloidal dispersions allows modifying the range 
and strength of attraction of the effective interactions between the colloidal particles. 
Adjusting the range of the attraction enables manipulating the topology of the phase 
diagram of a colloid-polymer mixture PQ 121 EH El- Both the nature of the demixed 
phases as well as the colloid and polymer concentrations at which demixing takes place 
depend on the range of attraction [S]. Industrially, it is relevant to understand the phase 
behaviour of colloid-polymer mixtures because colloidal particles and polymer chains are 
often jointly present in various products, such as food dispersions [HI Ej. Theoretically, 
the focus has been mainly on hard spheres dispersed in polymer solutions. This is due 
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to the fact that the hard-sphere dispersion is relatively simple and well-understood. 
Moreover, hard sphere-like colloids can actually be synthesised chemically. Examples 
are, for instance, stearyl-silica colloids |8 or PMMA-PHSA particles 

In practice, many stable dispersions containing spherical colloids consist of particles 
that are not 'pseudo-hard' but which can be characterised by a pair potential containing 
an additional soft repulsive tail. An example is a stable dispersion of charged colloids 
in an aqueous salt solution for which the interactions are described by the Derjaguin- 
Landau-Verwey-Overbeek (DLVO) theory ^Oj- This theory predicts that the effective 
pair interaction between charged colloids consists of a hard-core repulsion due to the 
finite size of the colloids and a screened- Coulomb (Yukawa) repulsion with the screening 
length given by the Debye length k~ x of the solvent. The screening length n~ l defines 
the thickness of the double layer of opposite charge surrounding each colloid. The range 
k _1 of the screened- Coulomb repulsion is a function of the salt concentration of the 
solvent, the dielectric constant and the temperature. Adjusting the salt concentration 
may influence the stability of a dispersion of charged colloids mixed with a neutral 
depletion agent in a common aqueous salt solution [HI HTJ [T3 [TH| . Hebert [TT] studied the 
precipitation of the charged, rod-like tobacco mosaic virus (TMV) by adding the neutral 
polymer polyethylene glycol (PEG). At similar PEG concentrations, precipitation of 
TMV was enhanced by adding salt. Patel and Russel [H] studied the phase behaviour of 
mixtures of charged polystyrene latex colloids and dextran as (neutral) polymer chains 
and reported a significant shift towards higher polymer concentrations of the fluid- 
fluid binodal curve as compared to predictions for neutral polymer chains mixed with 
hard spheres. Grinberg and Tolstoguzov jHj presented generalised phase diagrams of 
proteins mixed with neutral non-adsorbing polysaccharides in aqueous salt solutions. 
The miscibility or compatibility was shown to increase when the ionic strength of the 
solvent was lowered. The compatibility especially increased below 0.5 M. Finet and 
Tardieu ^2] studied the stability of solutions of the lens protein crystalline. Adding an 
excess of salt to this system does not destabilise the protein dispersion. Hence, it follows 
that the effective attractions between the proteins are absent or are very weak in the 
case of screened charges. Adding PEG however induces significant attractions ^2], and 
results in a shift of the liquid-liquid phase transition to higher temperatures ^Hj- Adding 
excess salt and PEG induces instant phase separation ^2|- A similar synergetic effect of 
salt and PEG was found in aqueous solutions of (spherical) brome mosaic virus particles 
In conclusion, the trend found in experimental studies on mixtures of charged 
'colloids' plus neutral polymers is that the miscibility is increased upon decreasing the 
salt concentration, i.e., increasing the range of the screened- Coulomb repulsion. 

In the light of these findings it is important to study theoretically mixtures of 
colloids with a screened-Coulomb repulsion mixed with neutral polymer chains and to 
investigate whether the trend found in many experimental studies is recovered. The 
amount of theoretical work performed so far is rather limited. Ferreira et al. [16 
made a PRISM analysis up to the level of the pair interaction and computed gas-liquid 
spinodal curves from the effective colloid-colloid structure factor. Denton and Schmidt 
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[T7] proposed a simple theory yielding the gas-liquid binodal curve. The fluid-solid 
coexistence curves were not considered and none of these theories were tested against 
computer simulations. Here we study the effect of a screened- Coulomb interaction on the 
total effective colloid-colloid interaction and on the resulting gas-liquid and fluid-solid 
phase transitions in a charged colloid dispersion with added non-adsorbing polymers. We 
demonstrate the fluid-solid coexistence is especially sensitive to the screened- Coulomb 
repulsion. The outline is as follows. First, in section |2] we calculate the total effective 
pair potential by juxtaposing the depletion and the DLVO (screened- Coulomb repulsion) 
contributions. In sectional a simple Ansatz for calculating the phase behaviour of the 
charged colloid dispersion is addressed, followed by an analysis of the effect of non- 
adsorbing polymer using free volume theory in section 14.11 The theoretical results are 
compared with Monte Carlo simulation results in section H~2l This communication will 
be summarised with some concluding remarks. 

2. The model 

We consider a suspension of charged colloidal spheres immersed together with non- 
adsorbing polymer in a common solvent. As the differences in length- and time-scales 
between the solvent molecules and the colloids and polymers are large, we can regard 
the solvent to be an inert continuum. Charged-stabilized colloidal suspensions can be 
described by a mixture of iV c colloidal particles with charge — Ze, with e the elementary 
charge, ZN C counterions with charge +e and by pairs of salt ions of charge ±e. It 
is convenient to consider the system in equilibrium with a salt reservoir with density 
p T s . However, computer simulations of this model are prohibited by slow equilibration, 
due to a large number of small counterions and salt pairs compared to the number of 
colloidal particles. A more coarse-grained model can be obtained by formally integrating 
out the degrees of freedom of the co- and counterions in the partition function ^H]. The 
effective pair interaction between the charged colloidal spheres with hard-core diameter 
<j c , reads 



where = |Rj — R 3 -| and Rj are the positions of the centres of the colloids. The range 



of the repulsive tail is set by the inverse Debye screening length ko c = ^S7ilBcr 2 p r s , 
which is related to the salt concentration p r s in the reservoir, and to the Bjerrum length 
Ib = e 2 /e s kT with k Boltzmann's constant, T the temperature, and e s the dielectric 
constant of the solvent. The strength of the repulsion of the repulsive Yukawa interaction 
is determined by the parameter e = [Zj{\ + Kcr c /2)) 2 l B /a c . We neglected the higher 
order, effective many-body terms which are insignificant in the high salt concentration 
regime that is considered in this work |18j . In this approximation, the effective potential 
(H) is the usual DLVO potential [TQl UHl 1201 IHJl 1221 1231 12^1 - In figure □ an example of 
a typical Yukawa repulsion is given by the dot-dashed curves for [3e = 20 and ko~ c = 100. 




) 



for < a c 

otherwise 



(1) 
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Here (3 = 1/kT. 

Different procedures [*""3 EH EE EH EOI exist that give the same functional form of 
the effective potential ((TJ), but with a density dependent Ra c = <J (4irl B a^(Zp c + 2p s )), 
where p c and p s are, respectively, the densities of colloidal particles and added salt ion 
pairs, in the system. The density p s depends on the colloidal density and is smaller 
then the reservoir density due to salt exclusion or Donnan equilibrium jSUESEH]- For 
the parameters used in this work, f3e = 20 and na c = 50, the effective k differ from 



the reservoir na c = JSttIbO'^pI by less then 1% for a colloid density near close-packing 
p s = 1.4<7~ 3 when one applies Donnan equilibrium. Hence salt partitioning is hardly 
perceptible for the salt conditions considered here. 

A dilute solution of polymers in a theta-solvent can be represented by (non- 
interacting) ideal polymers. In a theta-solvent, the attraction between the polymer 
segments compensates their excluded volume: the chains do not feel one another. The 
non-interacting polymer coils are excluded from the colloids to a centre-of-mass distance 
(cr c + cr p )/2 |""""f|. The chains can be interpreted as freely overlapping spheres jUHlEE]- 
The effective diameter of the polymer coil a v is twice the depletion thickness around the 
colloidal particle. The pair-wise potentials in this simple model read jHZ| 



Where Rj and are the positions of the centres of the colloids and the polymer coils, 
respectively, and r^ = [r< — r 3 -|. We mention that the so-called Asakura-Oosawa-Vrij 
(AOV) model which is a simple idealised model for colloidal hard spheres and 

ideal non-adsorbing polymer, is recovered by setting e to zero in the effective potential 



However, even within the context of this highly simplified model it is often more 
convenient to adopt a more coarse-grained view of the binary mixture by ignoring the 
degrees of freedom of the polymer coils and using polymer-mediated effective interactions 
between the charged spheres. Below we discuss the resulting effective interactions 
between the charged spheres due to the pair interactions of our model. 

By integrating out the degrees of freedom of the polymer coils, we can map 
the binary mixture of colloids and polymers onto an effective one-component system 
interacting with an effective one-component Hamiltonian "*"Hll39| 14*0*1 141| . This effective 
Hamiltonian consists of zero-body, one-body, two-body, and higher-body terms. In this 
work, we perform Monte-Carlo simulations for the effective Hamiltonian truncated at the 
pairwise term, and we determine the phase behaviour of the effective one-component 
system. We first derive an exact expression for the polymer-mediated effective pair 
potential or depletion potential. In general the effective pair potential can be calculated 
using the generalised Gibbs adsorption equation for two colloids separated by a distance 
Rij in a sea of polymer at fixed chemical potential p p , i.e., the system of colloids plus 
polymers is in thermodynamic equilibrium with a reservoir with only polymer chains 





for |Rj 
otherwise 




(2) 



(3) 
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(in a background solvent) at chemical potential \i v (321 EH1 113 EHj 

f3dn P {Y{R ij ,n P ) - Y(R ij = oo,/i p )) . (4) 

-oo 

The chemical potential is related to the fugacity of the polymer chains z p through 
Pdfi p = d\nz p . The quantity r(i?y,/i p ) is the adsorption or excess amount of polymer 
chains when the colloidal spheres are separated by a distance Rij. 

The depletion potential of the present model is identical to the depletion potential 
of the AOV model as the pair potentials U cp and U pp are the same in both models. Since 
the polymer chains do not interact, U pp = 0, we can replace z p by p p , where p p is the 
polymer density in the corresponding polymer reservoir. Further within the model the 
depletion layers have a fixed width of 5 = cr p /2 being close to the radius of gyration of 
a polymer chain EH Ej , so the excess adsorbed amounts are directly related to the 
overlap volume V over i ap (Rij) of depletion layers: 

I"(/t!jj. pjp) T(Rij OC. flp) P p Voverlap(Rij) 

A geometrical calculation for a c < R^ < a c + o~ p , the range where the depletion layers 
overlap, yields an analytic expression for V over i ap (Rij) and the depletion potential U^ ep 
reads for our model [SH ESI EE1 EI] 



PUfeJRij , Hp) 



™ 3 p p r p (1 + q) 3 



3 



3-R,;, R 



2(1 + q)a c 2(1 + qfa 3 c _ 
for cr c < R^ < cr c + a p (5) 
= for R^ > o c + Op 

where we define the size ratio q = o p ja c . This Asakura-Oosawa pair potential 
describes an effective attraction whose depth increases linearly with the polymer density 
in the corresponding reservoir p p . For convenience, we define the relative polymer 
concentration (jf p = 7Tcr p p p /6. Hence, <p p = 1 defines the overlap concentration of the 
polymer solution. The range of the potential is given by a p . The depletion interaction 
Udep(Rij) is plotted as the dashed curves for q = 0.1 and 0^=0.2 in figureHk, for q = 0.6 
and (fi p = 1 in figure HJo and for q — 1 and (fi p = 1.4 in figure [TJ:. By adjusting q and (f) p 
one can manipulate the range and strength of the depletion interaction. We chose 4> p 
such that PU^ e -p(cr c ) = — 3.5 for every value of q. 

It is important to note that for sufficiently large polymer coils, effective three- 
and higher-body interactions can not be neglected. More precisely, we expect an 
increasing number of higher-body interactions to become non-zero when q increases. 
It was shown |4Sj that the many-body character of the polymer-mediated effective 
interactions between the colloids yields a bulk phase diagram that differs substantially 
from those found for pair-wise simple fluids. For size ratios q < 0.1547, three- and many- 
body interactions are identical to zero and the mapping of the binary mixture onto the 
effective one-component Hamiltonian based on pair-wise additive effective potentials is 
exact. 
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Using the two-body approximation to the effective Hamiltonian, the total effective 
pair potential [/tot reads 

Utot(Rij) = U cc (Rij) + Ud ep (Rij)- (6) 

Examples of U^ t(Rij) are plotted in figure C] as the full curves, which are the sums 
of equations (JTJ) and (jSJ), denoted by the dot-dashed and dashed lines, respectively 
For q = 0.1 (figure^) there is a significant effect of the repulsive tail on the effective 
interaction as compared to the pure depletion contribution U^ ep (Rij). There is, however, 
still some attraction in the Ut t(Rij) curve between the charged repulsive spheres though 
it is significantly reduced as compared to the pure U^ ep (Rij) result. In figuresHJo and[Tt 
the main parts of the pair interaction curves Ut t(Rij) are identical to the pure depletion 
part U^ ep (Rij) for q = 0.6 and q — 1. Only the attraction at short interparticle distances 
is affected. The results plotted in figure [Hare qualitatively similar to the PRISM results 
of Ferreira et al. PB[ (see their figure 5). 



3. Phase behaviour of Charged colloids 

In this section we propose a simple description for the phase behaviour of charged 
colloidal suspensions. We consider the screened- Coulomb repulsion as a perturbation of 
the hard sphere interaction which is only valid for highly screened colloidal suspensions. 
The effective volume fraction 7] e of the charged repulsive spheres is approximated as 

Ve = (o^) Vc = mT]c ^ 

where r\ c = 7ia^p c /6 with p c the colloid number density and a e the effective diameter of 
the spheres defined by the Barker-Henderson relation [%9] : 

/•OO 

a e = a c + (1 - exp[-(3U cc (r)})dr (8) 

JtT c 

which is a useful way to generalise pair potentials of various shapes [SI], in particular, 
perturbations from hard-sphere behaviour (see [SID- Instead of U cc defined in equation 
one may use any other form for a (soft) repulsion. The physical effects are contained 
in the value for m. Since we assume that the collection of charged spheres behave 
similarly as a collection of pure hard spheres plus a small perturbation, we may use the 
Carnahan-Starling (CS) expression [H2] for the Helmholtz free energy to describe the 
thermodynamic properties of the fluid of charged spheres: 

mff uid (v c ,T) = Ve \n Ve + 4 ,f - ?| - Ve + Ve In ^§ (9) 

(1 - T] e ) 2 Vol 

where A c is the thermal wavelength of the colloids and where the hard sphere volume 
fraction r] c in the classical CS expression is replaced by the effective volume fraction 
rj e , defined in equations ([7J) and ©• In equation (JUJ) we use (as in [23]) the normalised 
Helmholtz free energy f c , defined as j3F c v c /V, where F c is the Helmholtz free energy, 
and v c = is the volume of a single colloid. We note that there are more advanced 



Phase behaviour of charged colloidal sphere dispersions with added polymer chains 7 



theoretical methods for Yukawa fluids j^lj that describe the simulation results for longer- 
ranged repulsions. 

The equation of state of the face-centred-cubic (fee) crystal phase of pure hard 
spheres is described accurately by the expression proposed by Hall j^S]. Likewise, the 
equation of state for the fee crystal phase of the charged spheres reads 



m/ crystal (r/c;T) = ^ 2 . 1306 + 3 ln 



Ve 



6A 3 

+ 77 e ln— f (10) 



.1 - Ve/Vcp 

containing the volume fraction at close packing i] cp = 7Ta/2/6 ~ 0.74. The quantity 
2.1306 is derived from computer simulation results [56] . 

Let us first check the accuracy of equations (jHI) and (fTTljl on the level of the resulting 
osmotic pressure Il c that follows from j3H c v c = r] c df c /dr] c — f c . The osmotic pressure 
for the fluid phase reads: 

m{3Tl^v c ( Vc ,T) = Ve + f + ^ 3 " ^ (11) 

(1 - Ve) 

and for the fee crystal, the pressure is: 

m(3 u c jy^v c ( Vc ,T) = _ 3r?e , (12) 

A Ve/ Vcp 

Results using equations (fTTj) and (|T2*j) with U cc (Rij) given by equation ((TJ) are plotted 
in figure El (curves) and are compared with computer simulation data (symbols), for 
(3e = 20 and 1/kct c equals (full curves, open circles), 0.01 (dotted curve, filled squares) 
and 0.02 (dashed curve, open triangles), corresponding to m = 1, 1.110 and 1.225, 
respectively. The results for pressures below (3H c v c = 6.2 correspond to the colloidal 
fluid phase, while the results for larger pressures correspond to the fee crystal. Figure El 
shows that the pressure increases upon increasing the range of the soft repulsion. The 
simulation results are well described by equation (fTTj) (see the results for /3Tl c v c < 6.2.) 
for the fluid phase for m — 1, m — 1.110 and m = 1.225. In addition, the results for the 
solid phase using equation (j!2)l ((3U c v c > 6.2) agree well with the simulations for m = 1, 
m = 1.110 and m = 1.225. Our results demonstrate that for sufficiently short-ranged 
soft repulsions, the screened- Coulomb interaction can be treated as a perturbation of 
the hard sphere potential using the Barker-Henderson relation (jHJ). 

The fluid-solid transition, first established for pure hard spheres by Alder and 
Wainwright and Wood and Jacobson [5E], can now be studied as a function of 
the softness of the repulsive tail. We determine the densities of the coexisting phases 
by equating the osmotic pressures and the chemical potentials = df c /drj c using 
equations and (fTU|). It becomes evident that the Ansatz of equations (J7J), ©, © 
and (fTTH) is expected to be useful only for short-ranged soft repulsions, so for o e ~ a c or 
l/na c — > 0. In order to test whether this approach is valid for relatively short-ranged soft 
repulsions, we compare the predicted fluid-solid transitions with computer simulation 
data for small values of l/na c . In figure El the fluid-solid binodals calculated for 
1/kOc < 0.04 are plotted as the full curves. In figure Ek, the (3e = 20 case is considered. 
The dots are computer simulation results by Hynninen and Dijkstra jSH]- Open symbols 
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are the pure hard sphere computer simulation results of Hoover and Ree [UUJ . For 
1/ko~ c < 0.02 the agreement is excellent though for larger Debye lengths deviations are 
found. This deviation is not surprising since the simple theoretical method is based on 
a perturbation from the hard sphere system and is, hence, expected to be only accurate 
for very small values of 1/ko~ c . Still, this simple theory for charged colloidal sphere 
dispersions suffices our purpose of studying the stability of charged colloid-polymer 
mixtures in the regime of l/no~ c < 0.02. To be more specific, this means that our 
Ansatz describes well the case of colloidal particles with a diameter of, say, a c = 100 
nm with a Debye length smaller than about 2 nm or, equivalently, an ionic strength 
> 0.02M (in the case of monovalent ions in water at room temperature). Actually, in 
nature, the ionic strength in aqueous salt concentrations is usually above 0.02M, so for 
large spheres we still capture the biologically relevant ionic strength regime. For globular 
proteins we can only compare with experimental data at significant salt concentrations 
but we can still make estimates of the main trends. 

In figure Eb, we compare our theoretical results (full curves) with computer 
simulation data (crosses) for (3e = 39. Again, our Ansatz is in fair agreement with 
the computer simulation results (data points) [SH]- Hence we validated our approach to 
describe the fluid and fee crystal equations of state and the fluid-solid transition of a 
charged sphere dispersion for 1/ko~ c < 0.02. This provides a base for studying the effect 
of adding non-adsorbing polymer to such a suspension. 



4. Phase behaviour of Charged colloid-polymer mixtures 

4-1. Free volume theory 

We study the stability of a mixture of charged colloidal spheres and non-adsorbing 
polymer chains in a common (background) solvent. A simple approach that successfully 
describes the stability of polymer-colloid mixtures is the semi-grand canonical free 
volume or osmotic equilibrium theory jHHSj. In this approach, a macroscopic volume V 
at temperature T is considered, that contains N c colloids together with polymer chains 
and solvent, which are in osmotic equilibrium with a reservoir containing only solvent 
and polymer chains. Hence, the system is considered in the (N c , V, z p , T) ensemble, in 
which the number of colloids N c and the fugacity of the polymer chains z p are fixed. 
Consider the thermodynamic identity 



-p 



,{d(3F(N c ,V,z p ') 



(3F(N C} V,z p ) = (3F(N c ,V,z p = 0) + ^ dz p ' ^ M ^ '■ PJ j (13) 

where we dropped the temperature dependence for convenience. The integrand can now 
be Taylor expanded about z p = 0: 

PF(N C , V, z p ) = PF(N C , V, z p = 0) + z p ( d ^ N ^ z A 

V P / z p =0 

+ 0(zl) + --- (14) 
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where the partial derivative can be written as 

(d/3F(N c , V, z p ) \ _(dF\ d/3fi p (N p ) Zp -_ 



V dz p J \dfi p J Zn=0 dz. 



z p 



(15) 



The number of polymer chains (N p ) Zp=Q can be related by definition to the averaged 
free volume that is available for the polymer chains (Vf ree ) Zp= o in the system of spheres 
that is undistorted by the addition of polymers: 

<i\y 2p=0 = P r p {v free ) Zp=0 . (i6) 

where p p is the density of ideal polymer in the corresponding reservoir. Defining the 
free volume fraction a = (Vf ree ) Zp =o/V , we can rewrite equation (JT3J) as: 

f3F(N c , V, z p ) = (3F(N C , V, z p = 0) - ap r p V + 0(z 2 p ) + ■■■ (17) 

The free volume theory [3] retains only the first-order term, neglecting terms 0(z p ) and 
higher. With this assumption the normalised thermodynamic potential / = (3Fv c /V 
can be written as the sum of two terms: 

f(N c , V, z p ) = f c (N c , V) - a4> r p q- 3 (18) 

The first term in equation ()18|) is the normalised Helmholtz free energy f(N c ,V,z p = 
0) = f c (N c , V) of a 'pure colloid' fluid at a given r] c , while the second can be interpreted 
as a perturbation due to the presence of polymer chains. Note that equation (fTH|) - (fTH|) 
holds for any colloid-polymer and polymer-polymer interactions within the assumptions 
that are made. However, for interacting polymers, a different approximation of the 
free volume theory gives a better description of the experimental phase diagrams [53J. 
All information about the interactions between colloid and polymer is contained in the 
variation of a with rj c . For our model and the AOV model, the free volume fraction 
a can be calculated accurately from scaled particle theory 0] (see e.g. Meijer jHS] 
for a comparison with computer simulation results). Once the coexisting colloid volume 
fractions are determined for given (jf pl the actual relative polymer concentrations can be 
obtained in the coexisting phases from <p p = a<p r p . An expression for the free volume 
fraction available for the polymer chains in a mixture of charged spheres using scaled 
particle theory reads 

a = (1 - Vc)exp(-A 1 - - C( - 3C( 2 - 3CC 3 ), (19) 

where 7 = 1] C /(1 — i] c ), ( = i] e /(l — i] e ), A = 3q + 3q 2 , B = 9q 2 /2, and C = q 3 . 
An explicit derivation is given in Appendix 1. Equation (fTT?|) reduces to the classical 
expression of Lekkerkerker et al. [Sj for the case r] e — r\ c (m = 1). The free volume 
available for the polymer chains is thus mainly a function of the pure hard sphere volume 
fraction r] c ] it is only affected by the screened- Coulomb repulsion at high colloid volume 
fractions or large polymer chains. In the derivation of equation (fT9*j) an equal statistical 
weight is assigned to all (non-overlapping) hard-sphere configurations, whereas the 
weight should involve the polymer-mediated effective interactions and the screened- 
Coulomb interactions. We expect our results to be accurate only if cr e /a c — 1 << q. If 
the depletion layers become small compared to a e — a c one expects hardly any overlap 
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of depletion layers. In figure 0] we compare the result of the free volume fraction a as a 
function of the colloid volume fraction r] c of equation (fT9"|) with Monte Carlo simulation 
results for (p P ) Zp / p r p using equation (|%2*|) for q = 0.1. It is worth mentioning that a is 
evaluated in the pure charged colloid system in the free volume theory, i.e., z p = 0, 
while (pp)z p /Pp from simulations do depend on z p . In the simulations, however, we 
use z p along the bulk binodals as shown in figure [Hi (which will be discussed later). 
The technical details of the simulations are described in Appendix 2. The theoretical 
curves in figure 0] show only a slight effect of the screened- Coulomb repulsion on a for 
r] c > 0.4. Within the statistical error bars, no significant effect of the screened-Coulomb 
repulsion was detected in the simulations on the free volume fraction. Moreover, the 
^-dependence is hardly noticeable. We thus conclude that equation (fT9"]) is accurate for 
short-ranged screened-Coulomb repulsions for q = 0.1 or, equivalently, that the effect of 
the screened-Coulomb repulsion on a is negligible. 

We can now analyse the effect of the screened-Coulomb repulsion on the phase 
behaviour. The osmotic pressures II and colloid chemical potentials p, c can be found by 
differentiation of the Helmholtz free energy ()18|) . The colloid volume fractions in each 
of the coexisting phases, r}\ and rj^, are obtained by equating \i c and II at fixed polymer 
reservoir concentration (ff . The Helmholtz free energy density (jlOj) is used for the solid 
phase, while equation @ is employed for the fluid. The free volume theory predicts for 
the AOV model, i.e., a system of colloidal hard spheres and ideal polymer, a broadening 
of the fluid-solid transition with increasing <jf for size ratios q < 0.3, while the fluid-fluid 
transition is metastable with respect to a broad fluid-solid transition |3| I41j. 

We now turn to our case of charged colloidal spheres and ideal polymer chains. In 
figure we plot the phase diagram using the free volume theory (|T8|) for a size ratio 
q = 0.1 and (kcxJ- 1 = 0.0, 0.005, 0.01, and 0.0125 in the (ij e , (j> r p ) plane. For {na c y l = 
and (ff p = 0, we recover the well-known pure hard sphere freezing transition at r\ c = 0.494 
and 0.545 |0Dl- in the case of charged spheres, the freezing transition at (j) r p = shifts to 
lower colloid volume fractions r] c , which is in line with the results in figure El and is due 
to a larger effective volume of the charged spheres. Figure El shows clearly that the fluid- 
solid transition widens upon increasing the polymer concentration. More specifically, 
the broadening of the freezing transition shifts to higher 0£ with increasing range of the 
screened-Coulomb repulsion (/tcr c ) _1 . This can be explained as follows. Upon increasing 
the range of the screened-Coulomb repulsion, i] e = mrj c , and hence f c , increases. At the 
same time the free volume fraction a is not affected significantly (see figure HJ upon 
adding a screened-Coulomb repulsion. So, in order to attain a similar effect on / (see 
Eq. (18)), a higher polymer concentration is required to broaden the freezing transition. 

For larger values of q, say q > 0.4, a fluid-fluid coexistence becomes stable in 
the AOV model, which dominates the phase behaviour at colloid volume fractions 
rj c < 0.49. In analogy with figures [Ud and Ht we choose q = 0.6 and 1 and study 
the effect of the repulsive screened-Coulomb interaction on the fluid-solid and fluid-fluid 
transition. In figures El and we plot the predictions from free volume theory in the 
(rj c , <f>p) representation for q = 0.6 and q = 1.0, respectively. We again find that the 
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freezing transition at <fi r p = shifts to lower colloid volume fractions r] c upon increasing 
(kct c ) -1 . In addition, figures |H] and UJ show that the fluid-fluid demixing shifts to higher 
with increasing range of the screened- Coulomb repulsion (kct c ) _1 . Hence, the trends are 
similar as for q = 0.1. The screened- Coulomb repulsion reduces the depletion effect. The 
critical points are indicated as the filled circles in figure EH and and they indicate that 
the critical colloid volume fraction r] c shifts to somewhat smaller values upon increasing 
the range of the soft repulsion. 

Finally, we convert the polymer reservoir concentration 0£ to that in the actual 
system p . Figure |Hl shows the conversion of the phase diagram of figure El for q = 0.1 
and f3e = 20 and Debye screening lengths (ko" c ) _1 = (dashed curves; the pure hard 
sphere case), (kUc)" 1 = 0.005 (dot-dashed curves), (k(j c ) _1 = 0.01 (full curves) and 
0.0125 (dotted curves). The phase stability of a mixture of charged colloids and neutral 
polymer chains in an aqueous salt solution is thus expected to depend very sensitively 
on the screening length, and thus on the salt concentration, at least for small size ratios 
q. In figure EH we investigate the effect of (3e on the phase behaviour. We plot the 
converted phase diagram for the same set of parameters as in figure |H1 i.e., q — 0.1 and 
varying Debye screening lengths (kct c ) _1 , but with f3t = 39 instead of (3e = 20. There 
is a striking similarity between the two sets of results and, hence, we conclude that the 
effect of (3e is not significant. For larger values of j3e, the system becomes more sensitive 
to ko c and thus to the salt concentration. Finally, figure ITUl shows the conversion of the 
phase diagram as shown in figure UJfor q = 1.0 and f3e = 20 and Debye screening lengths 
(kct c ) _1 = (dashed curves; the pure hard sphere case), (na c )~ l = 0.01 (full curves) 
and 0.02 (dotted curves). 

It follows that an increase of the reduced Debye length (kct c ) -1 shifts the fluid-fluid 
coexistence curves upwards. Using a PRISM approach, Ferreira et al. [T0J also found 
this trend for the spinodal curve of the demixing fluid (see their figure ITTj) based on 
determining the composition where the inverse structure factor vanishes in the long 
wavelength limit. Figures 171 and fTUl show clearly that the shift in polymer concentration 
of the fluid-fluid binodals for q — 1 is weak compared to the shift in the fluid-solid 
binodals for small q. 

A relevant quantity that measures the relative influence of the screened-Coulomb 
repulsive pair interaction is (/tcr p ) _1 or (Kcr c g) _1 . Hence, the size of the polymer chains 
(or the depletion thickness) as compared to the range of the repulsion determines 
the relative importance of the screened-Coulomb repulsion on the total effective 
depletion interaction. In biological systems such as charged proteins mixed with 
neutral polysaccharides, where often k~ x < a c q, we expect that the phase behaviour 
is only moderately sensitive to the salt concentration. Decreasing the salt concentration 
significantly is then expected to stabilise the charged biocolloid dispersion against 
depletion-induced demixing. This explains the enhanced miscibility found in mixtures 
of proteins mixed with neutral non-adsorbing polysaccharides in aqueous salt solutions 

In many charged colloidal dispersions the soft repulsion is expected to suppress 
the depletion effect. In several applications such as paints and food dispersions where 
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colloidal particles are mixed with polymer chains, a screened-Coulomb repulsion helps 
stabilising the dispersion. 

4-2. Simulations using the effective Hamiltonian 

In order to test the validity of the predictions from free volume theory, we compare our 
results with Monte Carlo simulations. We determine the phase diagram of the effective 
one-component system by calculating the dimensionless free energy density / = (3Fv c /V 
as a function of the colloid packing fraction rj c and the fugacity of the polymer chains z p 
in simulations. For non-interacting chains the fugacity z p equals the density of polymer 
chains p r p in the corresponding reservoir. As the free energy density cannot be measured 
directly in a Monte Carlo simulation, we use thermodynamic integration to relate the 
free energy of the effective system to that of a reference system at the same colloid 
volume fraction rj c . To this end, we write the total free energy density as the sum of 
two contributions 



where f c (N c , V) is the free energy density for a system of N c charged colloids in a 
volume V interacting with a hard-core repulsive Yukawa potential ((TJ), and /^ep * s the 
contribution of the depletion potential (jlj]) to the free energy density. The free energy 
density f c {N c , V) is computed using the so-called A-integration for the fluid phase and 
the standard Frenkel-Ladd method for the solid phase [HI] with the Einstein crystal 
as a reference. A second A-integration is then carried out, for both the fluid and the 
solid phase, to determine the free energy density contribution /d ep . For more details 
of the simulations, we refer the reader to Appendix 2. In order to map out the phase 
diagram we determine the total free energy density f(r] c , z p ) for many state points (r] c , z p ) 
in simulations. We employ common tangent constructions at fixed z p to obtain the 
coexisting phases [10] ■ In order to assess the performance of the free volume theory, 
phase diagrams were determined for some of the parameters as in figure El (q = 0.1). 
The results are shown in figure ^2 and agree semi-quantitatively with those in figure 
The main difference is due to the fact that the results of both approaches deviate 
already for the AOV model (hard sphere colloids with ideal polymer, i.e., (k<t c ) _1 = 0) 
[¥T] . For the AOV model, the theoretical binodal is shifted with a factor of about 1.3 
in (f)^ compared to the simulations. This factor between the theoretical predictions and 
the simulation results is about 1.2 when the screened-Coulomb repulsion is added. 

We also compare the free volume theory results for the phase behaviour for q = 0.6 
and 1 as shown in figures |H1 and [7| with computer simulations. The phase diagrams 
obtained from Monte Carlo simulations of the effective one-component systems are 
plotted in figure (q — 0.6) and figure (g = 1.0). The main theoretical trends 
are also found in our simulation results. Again, the data of the simulations suggest 
that a higher polymer concentration is required to induce the fluid-fluid transition upon 
increasing the range of the soft repulsion. 




(20) 
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Finally, we stress that the free volume theory incorporates some of the many-body 
effects which are present at large q, while our simulations are based on a two-body 
approximation to the effective Hamiltonian. It is therefore difficult to make a direct 
comparison between the simulation results and those obtained from free volume theory. 
However, for q < 0.1547, the mapping of the charged colloid-polymer mixture onto an 
effective one-component Hamiltonian based on depletion pair potentials is exact and 
thus a direct comparison is feasible for our results for q = 0.1. 

5. Concluding remarks 

We have studied the effect of a short-ranged screened-Coulomb repulsion on the phase 
stability of mixtures containing charged spheres and non-adsorbing polymer chains. The 
charged spheres are described as hard spheres with an additional screened-Coulomb or 
Yukawa repulsion with the screening length given by the Debye length setting the 
range of the soft repulsion. The phase behaviour of the charged sphere dispersion is 
described using standard expressions for the colloidal hard sphere fluid and fee crystal 
with the hard sphere volume fraction replaced by an effective volume fraction that de- 
pends on the Yukawa interaction between the spheres. Our results obtained from free 
volume theory and Monte Carlo simulations show that the additional screened-Coulomb 
repulsion reduces the depletion effect. For mixtures of small polymers plus relatively 
large charged spheres the fluid-solid transition is shifted to significantly larger polymer 
concentrations with increasing Debye screening length while for relatively larger 
polymers the effect is weaker: the resulting fluid-fluid binodal is affected weakly by 
adding a short-ranged soft repulsion. In general, the range of the screened-Coulomb 
repulsion compared to the range of the depletion attraction determines qualitatively the 
reduction of the depletion effect, and hence, shifts the fluid-fluid and fluid-solid binodals 
correspondingly. In conclusion, a mixture of charge-stabilised colloids and non-adsorbing 
polymers at large concentrations of both components can be stabilised by lowering the 
salt concentration, which increases the range of the screened-Coulomb interaction of the 
colloids. 
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Appendix 1: Free volume fraction in a charged sphere dispersion 

In this appendix we consider the free volume fraction a that is available for ideal polymer 
chains in a sea of charged spheres with diameter o c . As the centre-of-mass of the 
polymer chains is excluded from the centre-of-mass of the charged colloids by a distance 
(<r c + o~ p )/2 and the polymer interactions are ideal, a is just the free volume fraction 
for a single hard sphere with diameter a p in a sea of charged spheres. This free volume 
fraction a can be determined from the chemical potential of the polymer chains. The 
chemical potential for inserting a polymer in a sea of charged spheres consists of an ideal 
gas term and a work term W. 

Pfx p = \np p Al + W (21) 

Following Widom's particle insertion method [65J, the required work to insert a polymer 
W is equal to (3W = — In a. The work W can be determined from scaled particle 
theory, that considers two limits. In this theory, the size of the particle is scaled with 
a parameter x. In the limit x <C 1, the polymer coils reduce to points, and hence the 
volume fraction available to the polymer is simply unity minus the sum of the overall 
hard sphere volumes plus the depletion layers around them: 

a{x < 1) = 1 - ^ Pc (a c + xa p f (22) 
o 

where p c = N c /V is the number of colloidal spheres per volume (related to rj c via 
rj c = np c al/Q). Hence, it follows 



(3W{x < 1) = -In 



7T . .3 
-Pc(0"c + XCTp) 
D 



(23) 



On the other hand, if x ^> 1, the work required to insert a polymer coil in a sea of 
charged spheres, is approximately the work to create a hole with the size of the polymer 
coil, which is equal to the volume of the polymer coil times the osmotic pressure II C of 
the dispersion of charged colloids: 

W(x > 1) = ^x 3 a 3 Jl c . (24) 



,.2 



In scaled particle theory, W{x <C 1) is expanded about x = up to order x and 
W(x ^> 1) is added as the x 3 term. 

w(x) = w{x = 0) + x(f + i^(^)^ + ^(25) 
Scaling the polymer coils to the desired size by x = 1, yields 

(3W{x = 1) = - In [1 - r) c ] + 3g 7 + - (V 7 + 9g 2 7 2 ) 

+ \^ p (3Yi c (26) 

where 7 = r) c /(l — r] c ). Hence a follows straightforwardly from W = W(x = 1). For 
pure hard spheres one usually takes the Percus-Yevick result for the pressure Hhs from 
the virial route (see [HI]) since it is consistent with SPT 

P^hs 1 + Vc + Vc 1 3t? c 3 V 2 c 

Pc (1 - Vcf 1 - Vc (1 - Vcf (1 - Vc)' 
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Inserting this expression for Hhs into equation |2H1 for Il c yields 

/3w = - In [1 - Tic] + (A + C) 7 + (B + 3C) 7 2 + 3C 7 3 (28) 

where A, B, and C, are defined below equation Hence, we arrive at the standard 
SPT result for the free volume fraction of ideal polymer in a sea of hard spheres [Sj 

[(A + C) 7 + (B + 3C) 7 2 + 3C 7 3 )] (29) 

In case of colloidal spheres interacting with a Yukawa pair potential, we rewrite 
equation (j27j) following the approach outlined in section 3 giving a pressure 

P^Yuk = 1 + Ve + V 2 e = 1 % 3^ 

Pc " (1 - ^ e ) 3 " 1 - Ve (1 - r/ e ) 2 (1 - Ve ) 3 1 J 
where r] e is defined by equation (J7J). Using this expression for Il c in equation l2*o1 yields 
equation [THl 



a = (1 -?7 c Jexp 



Appendix 2: Technical details of the simulations 

This section describes the technical details of the simulations. To this end we consider 
the total effective one-component Hamiltonian of the colloids at fixed polymer fugacity 

Zp 

H(z p ) = H cc + H dep (z p ) (31) 

where the Hamiltonian H cc consists of a sum of colloid-colloid pair potentials U cc , which 
can be split into a sum of hard-sphere potentials Uhs an d a sum of repulsive Yukawa 
potentials: 

N c N c N c 

Hcc = ^2 Ucc{Rij) = ^2 ^Hs{Rij) + ^2 ^Yuk{Rij) (32) 

i<j i<j i<j 

with 

it it} \ exp(Ka c ( J R ij /a c - 1)) 

U Yu k{Rij) = e — -. for Rij > a c (33) 

Kij/a c 

and -ffdep a sum °f depletion potentials t/^ep © : 

N c 

H dep(zp) = J2 U dep(Rif,z P )- (34) 

i<j 

The different Hamiltonians give rise to the corresponding free energy contributions of 
equation (I2U1) . 

To determine the free energy contribution of the Yukawa potential, for the fluid 
phase, we introduce the auxiliary Hamiltonian 

N c Nc 

<f = E V H siM + A E U Yv k(Rij) (35) 

i<j i<j 

where < A < 1 is a dimensionless coupling parameter: at A = the auxiliary 
Hamiltonian is that of a pure system of N c hard spheres, while at A = 1 it is the 
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Hamiltonian of N c charged spheres. The free energy is determined by applying the 
standard A-integration technique [HI] 

f c (N c ,V) = f HS (N c ,V,\ = 0) + ^ f'dx/^pUYukiRij)) ■ (36) 

\*<i 1 n c ,v,\ 

The angular brackets denote a canonical average over a system of N c particles interacting 
with the Hamiltonian H^ d , while fiis(N c , V, A = 0) is the free energy of a system of 
hard spheres, for which we use the Carnahan-Starling expression [52J. We start the 
canonical simulations from a random, non-overlapping configuration, and use 15000 
MC sweeps per particle for equilibration and typically 15000 production moves for each 
value of the coupling parameter A. In principle, the free energy of the solid phase can 
be computed with the same technique using the Hall expression for the free energy of 
the hard sphere crystal [22] ■ However the latter is only properly defined for packing 
fractions larger than the value at hard sphere freezing i] c = 0.545, while charged spheres 
can yield crystal phases at lower packing fractions. A different technique is used for the 
solid phase by introducing the auxiliary Hamiltonian 

Nc N c 

H S S = EM%) + \k B T^-r ,i) 2 /o- 2 c , (37) 

i<j i=l 

where r o i denote the ideal lattice position of particle % in a fee crystal. The free energy is 
computed by applying the integration technique introduced by Frenkel and co-workers 



fc(N c ,V) = f2 4 (N c ,V,X = X max ) + f corr (N c ,V) 



i ein 

CM 



- i l dX (5> - r 0)i ) 2 /- c 2 ) , (38) 

where the angular brackets denote a canonical average of the mean square displacement 
of N c particles interacting with the Hamiltonian H^fl d , while the superscript CM 
denotes that it is calculated for a crystal with fixed centre of mass. The parameter X max 
is chosen such that for A = X max the system behaves like a non-interacting Einstein 
crystal with fixed centre of mass and Madelung energy t/y«fe(i"o °); i- e -5 the potential 
energy of a crystal with all particles at their ideal lattice positions. Typical values for 
Xmax range from 1000 to 100000 for high densities. The free energy of a non-interacting 
Einstein crystal with fixed centre of mass reads 



f™(N c , V, X = X max ) = pU Yuk (r»<) - ln 



7T 



A, 



(N c - l)v c , 
+ V — ' In 



A 3 . 



(39) 



The correction term f corr arises when the constraint on the centre of masses is released, 
i.e., the Helmholtz free energy difference between the unconstrained and constrained 
crystal: 



f corr (N c ,V) = ^\n 



A 3 



VNl' 2 



(40) 
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The equilibration is done for 15000 MC steps per particle, and the averages are taken 
for 15000 MC steps per particle. 

To determine the free energy contribution of the AOV depletion potential f^ep we 
employ a second thermodynamic integration for the solid and fluid phase: 

fde P = f(N c ,V, P ;)-f(N c ,V,p r p = 0) 

J0 V UP P > N c ,V,pr> 

The system at p r p = is a system of colloids interacting with pair potentials U cc . The 
integral is calculated by dividing the interval [0, p r p ] in 20 to 30 equally spaced intervals. 
We used 10000 MC steps per particles for equilibration, while averages were taken for 
20000 MC steps per particle. In addition we used the integrand to determine the number 
density of ideal polymer in the system using the thermodynamic relation 

/ df(N e ,V,fQ \ = f df(N c ,V,z p ) \ = _ v (_p^ 

V d P r P ) \ dz P ) Vc p r P 
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Figure 1. Pair interaction between two colloidal spheres with hard-core diameter <j c 
interacting with a Yukawa pair potential U cc (equation ^) and an effective depletion 
interaction t^jep (equation^ due to the presence of non-adsorbing polymer coils with 
an effective diameter a p . The total effective interaction U^ ^ (equation |SJ is denoted 
by the full curve. The Yukawa repulsion is characterised by /3e = 20 and na c = 100. 
The depletion interaction are: a) for a size ratio q = a p /a c = 0.1 and relative polymer 
concentration 0^=0.2, b) q = 0.6 and <£=1, and c) q = 1 and </ff=1.4. 
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Figure 2. Pressure of a dispersion of charged colloidal spheres interacting with 
U cc (Rij) ( see equation for (3e = 20 and 1/kct c = (open circles, full curves), 0.01 
(filled squares, dotted curves), and 0.02 (open triangles, dashed curves), corresponding 
to m = 1,1.110, and 1.225, respectively, in equation JJJ). The symbols denote the 
simulation results, while the curves denote the theoretical predictions of equations 
(ITTl (lower set of pressures) and fT2"|l (upper set). 
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Figure 3. Fluid-solid (fee) transition of charged colloidal spheres interacting with 
a hard-core repulsive Yukawa potential JU with a) (3e = 20 and b) f3e = 39. Open 
symbols denote simulation results of pure hard spheres taken from Ref. [HO] ■ Filled 
circles denote the simulation results of Hynninen and Dijkstra Full curves 

correspond to the theoretical predictions as described in section |3] 
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Figure 4. Free volume fraction a = (Vf ree ) z —o/V for a mixture of charged colloids 
and ideal polymer with size ratio q — a p /a c — 0.1 as a function of hard-core 
volume fraction r/ c . The screened- Coulomb repulsion is characterised by /3e = 20 
and various values of k<j c . Full, dotted and dashed curves represent equation H19(l for 
kct c = oo (m = 1), 100 (m = 1.110) and 80 (m = 1.138), respectively. Note that the 
differences in the theoretical curves are only noticeable at high r/ c . The symbols with 
errorbars denote the Monte Carlo simulation results for (p P ) Zp /Pp using equation [fl^ 
for kct c = oo (closed squares), na c = 100 (open circles) and na c = 80(closed triangles), 
where we used z p at bulk coexistence (see figure ITTJ) . 
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Figure 5. Phase diagram of a mixture of charged spheres and ideal polymer as 
obtained from free volume theory for q = 0.1 as a function of the colloid volume 
fractions t] c and the ideal polymer reservoir concentration <jf p . The screened-Coulomb 
repulsion (JTJ is characterised by j3e = 20 and various values of kct c as indicated. F 
and S denote the stable fluid and solid fee phase. F + S denotes the stable fluid-solid 
coexistence region. 




Figure 6. Same as figure but for q = 0.6. F + F denotes the stable fluid-fluid 
coexistence region. 
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Figure 7. Same as figure El but for q = 1.0. F + F denotes the stable fluid-fluid 
coexistence region. 
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Figure 10. Same as in figure |H1 but with size ratio q — 1.0. The filled circles indicate 
the location of the gas- liquid critical points. The inset is a blow-up of the critical 
region of the liquid-gas binodal. 
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Figure 11. Phase diagram of a mixture of charged spheres and ideal polymer as 
obtained from simulations of the effective pair potential J^J for g = 0.1 as a function 
of the colloid volume fractions r\ c and the ideal polymer reservoir concentration <$> r . 
The screened-Coulomb repulsion QJ is characterised by j3e = 20 and various values of 
na c as indicated. The curves serve as a guide to the eye. Open circles correspond to 
the pure hard sphere case (kct c ) _1 = 0, closed circles to a screened-Coulomb repulsion 
with (kct c ) _1 = 0.01 and open triangles to (kct c ) _1 = 0.0125. F and S denote the 
stable fluid and solid fee phase. F + S denotes the stable fluid-solid coexistence region. 
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Figure 12. Same as figure ITTI but for q = 0.6. F + F denotes the stable fluid-fluid 
coexistence region. 




Figure 13. Same as figure ITTI but for q = 1. F + F denotes the stable fluid-fluid 
coexistence region. 



